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Abstract 

We investigate the possible regular solutions of the boundary Yang-Baxter equation for 
the fundamental U q [G<2\ vertex model. We find four distinct classes of reflection matrices 
such that half of them are diagonal while the other half are non-diagonal. The latter 
are parameterized by two continuous parameters but only one solution has all entries 
non-null. The non-diagonal solutions do not reduce to diagonal ones at any special limit 
of the free-parameters. 
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1 Introduction 



It is, by now, well known that integrable two-dimensional systems of statistical mechanics can 
be derived from the Yang-Baxter equation [1]. The respective Boltzmann weights can be related 
in a natural manner to the elements of a Yang-Baxter R- matrix solution R a b(x) £ <8> 
satisfying the relation, 

Ri2(x)R 23 (xy)R 12 (y) = R 23 {y)Ru{xy)R 23 {x) (1) 

for arbitrary spectral parameters x and y. 

For an open statistical system not all possible types of boundaries are compatible with the 
Yang-Baxter condition (1). However, the integrability at the boundary may be assured when 
the boundary weights K(x) G C N fulfill the so-called reflection equation [2, 3], which reads 

R 12 (x/y)K(x)R 12 (xy)K(y) = K(y)R 12 (xy)K(x)R 12 (x/y) (2) 

i 

where K(x) = K(x) <g> In and In is the N x N identity matrix. 

An important class of _R-matrices are those based on a quantum group U q [Q] g-deformation 
of a classical Lie algebra Q [4]. For these models, from a given solution K(x) of Eq.(2) one 
can, in principle, construct families of commuting transfer-matrix [3, 5]. This fact has motived 
several authors, see for instance refs. [6, 7, 8], to search for X-matrices solutions associated to 
such quantum integrable models as well as for their elliptic extensions [9, 10]. More recently, 
attempts to classify all the ^-matrices of the non-exceptional U q [Q] vertex models either by 
direct analysis of Eq.(2) [11] or by others potentially systematic approaches [12, 13] have been 
discussed in the literature. 

In spite of all these works not much is known about the structure of the i^-matrices of the 
vertex models based on the quantum exceptional Lie algebras. In fact, even for the simplest 
such case, i.e. the fundamental U q [G 2 ] _R-matrix [14], only the diagonal solutions have been 
studied [15]. The purpose of this paper is to start to bridge this gap by presenting the complete 
reflection matrices associated to the minimal U q [G 2 ] vertex model. 
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This paper is organized as follows. In next section we describe the seven-dimensional U q [G?\ 
.R-matrix on the Weyl basis. This step makes it possible to adapt the method developed for 
non-exceptional models [11] to deal with the exceptional U q [G?\ case. In section 3 we discuss 
what we hope to be the complete set of reflection matrices. We find two families of non- 
diagonal solutions and we confirm as well the special diagonal ones given before in ref. [15]. We 
observe that the two types of diagonal .fT-matrices cannot be obtained as special limits of the 
non-diagonal ones. Section 4 is reserved for our conclusions. In Appendix A we present the 
explicit Boltzmann weights expressions of the U q [G 2 ] vertex model. 

2 The U q [G 2 ] R- matrix 

In this section we shall present the U q [G2\ -R-matrix in a suitable basis for the analysis of 
the reflection equation (2). In terms of the quantum group framework [4], this matrix can be 
written as linear combination of the U q [G?\ projectors operators Py A [14]: 

E Pa(x)P Va (3) 

A=0,Ai,A 2 ,2A 2 

where V\ denotes the irreducible representations occurring in the seven-dimensional U q [G?\ 
decomposition V\ 2 <g> V\ 2 = V ® V Al © V\ 2 © V 2 \ 2 . 

The weights p\(x) are functions of the quadratic Casimir element in an irrep with the 
highest weight. They satisfy the following relations, 

[21 [121 

PaM = 1 ix Po(x), P2A 2 (x) = [2} x [12} xPo (x), Pa 1 (x) = [12} xPo (x) (4) 

where [a] x = X ^,Z X an d Po( x ) is an arbitrary normalization. 

A direct analysis of the reflection equation (2) is in principle easily made by representing the 
K- matrix K(x) in terms of the Weyl basis. In this approach, the same has to be done for the 
U q [G2\ -R-matrix as well. This last step involves a considerable amount of additional work even 
when the projectors Py A are known in terms of the so-called g-Wigner coefficients. Ommiting 
here such technicalities, we find that the _R-matrix defined by Eqs.(3,4) can be rewritten in 
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terms of the following expression, 



R(x) 




a,/3 = l,^4 



<X<P,/3' 



7 



+ X! ^2(0;) [e aa ® e 44 + e 44 ® e a / a /] + ^(a) [e 44 ® e aa + B a > a > 59 e 44 | 



a=l 



o<4 
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+ H C(a,P) [e a/ 3 ® ep a + ep a <g) e a p + e a 'p> <g) ep> a > +ep>a> <8) e a ^] 



a, 0=1,^4 



a</3,/3> 
3 



+ Y {%] [ e «/3 ® e4<5 + e /3" ® e<54 + e<5 ' 4 ® e P' a ' + e4<5 ' ® e "'/ 3 '] 



fc=i r fc 



+ ^4[fe] [e Q '/3' ® e 4( y + e^/ ® e^ 4 + e,5 4 <g> + e 4(5 ® e a , 



(5) 



where a' = 8 — a and the symbol denotes the sum over a set of three indices {(a, (3,5)} 



such that r\ = {(1, 3, 2), (5, 2, 7), (6, 3, 7)}, T 2 = (1, 2, 3) and T 3 = {(2, 1, 5), (3, 1,6)}. We also 



recall that e a ^ refers to the standard 7x7 Weyl matrices. 

The explicit expressions of all the weights appearing in expression (5) have been summarized 
in Appendix A. An advantage of the representation (5) is that it exhibits explicitly the two 
£7(1) symmetries of the G 2 Lie algebra. This property will be very useful in next section in 
order to classify the independent functional relations for the reflection matrices. We also note 
that the last term in Eq.(5) represents additional Boltzmann weights as compared with those 
present in the i?(x)-matrix of the non-exceptional vertex models [4]. 

We close this section by mentioning useful relations satisfied by the matrix R(x) = PR(x), 
where P is the seven- dimensional permutator. Besides the standard properties of regularity 
and unitarity the matrix R(x) satisfies the so-called PT-symmetry given by 




(6) 



where the symbol tk denotes the transposition in the space with index k. 
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Yet another property is the crossing symmetry, 

Ri2(x) = ^Ki&fo 12 /*)^- 1 (7) 
where V is an anti-diagonal matrix whose non-null elements are Vn = 1,V26 = — V35 = 



3 The U q [G 2 ] K- matrices 



The purpose here is to search for the complete set of regular fT-matrices for the i/gfGy vertex 
model defined in the previous section. More specifically , we are interested in reflexion matrices 
having the general form, 

7 

K(x) = k a ^(x)e a>p (8) 

a,P=l 

with the constraint k a ^{l) = 5 a ^. 

Direct substitution of the ansatz (8) and the i?-matrix (5) in Eq.(2) give us several in- 
dependent functional equations for the elements k a ^(x). In order to select and solve these 
equations we act on them with the operator and afterwards we take the y = 1 regular limit. 
Here we shall denote these resulting algebraic equations by the symbol E[i,j] where the index 
refer to the ith row and jth column of the original reflection equation (2). Note that such 
equations involve only the variable x and a number of free-parameters ui a ^ defined as 

Uaj, - —jjj-\y=l (9) 

Our next step is to organize the algebraic equations E[i,j] in suitable blocks of relations 
B[i,j] involving related ii'-matrices elements. We find that such blocks are made by combining 
for a given pair [i,j] the E[i,j], E[j, i], E[50 — i, 50 — j] and £"[50 — j, 50 — i\ set of relations. 
Considering that we are looking for general non-diagonal ii'-matrices we start our analysis by 
inspecting the blocks B[i,j] possessing only off-diagonal elements k a ${x). The simplest ones 
are £[1,33], £[1,41], £[2,42], £[3,35], £[9,17] and £[12,20], providing us constraints to the 



elements k a>a i(x), which reads 

k a ,a'(x) = ^-k h7 (x) (10) 

^1,7 

By the same token, the blocks _B[1,47], S[l,48] and 5 [8, 38] lead us to determine the 
following elements, 



ki, 5 (x) 


= r(x) 


(ai (x) en (x) wi >5 


/ \ / \ \ ^1,7 \ x ) 

C 2 fx) C13 fx) W 3>7 ) 

^1,7 


(11) 


kz, i fa;) 




ai fx) en fx) ua 1 


— C? fx) Cis fx) Cc>7 q) 7,1 ^ ^ 
z w io v / ')*-*/ 

^7,1 


(12) 




= T(x)( 


61 (x) Cn (x) W 7j3 


/ \ / \ \ ^7,1 (#) 

c 2 fx) C31 fx) CJ 5j i) 

^7,1 


(13) 


k 3 ,7(x) 


= T(x)( 


6l (x) Cn (x) CJ 3i7 - 


C2 [X) C31 (x)cJi i5 ) 

^1,7 


(14) 


k 16 (x) 


= r(x) 


(ai (x) cn (x) wi >6 


/ \ / \ \ ^1,7 (^J 

c 2 fx) C12 fx) ^2,7) 

^1,7 


(15) 


kei(x) 


= r(x) 


[ai (x) cn (x) a; 6j i 


/ \ / \ \ ^7,1 0*0 

- C 2 fx) C12 fx) U 7 ,2) 

^7,1 


(16) 


k^ix) 


= r(x) 


'b^ fx) Ci 1 fx) too 7 


/ \ / \ \ki,7 \ x ) 

— Co I X ) C01 I X ) (Vi k ) — : 

^1,7 


(17) 


k 72 (x) 


= r(x) 


(61 (x) cn (x) uj 7i2 


/ \ / \ \ ^7,1 (#) 
C 2 fx) C21 f^)^6,l) 

^7,1 


(18) 


k 23 (x) 


= r(x) 


(ai (x) Cn (x) cj 2j3 


/ \ / \ \ ^1,7 \ x ) 

c 2 fx) c 25 (aOw 5)6 ) 

^1,7 


(19) 


h 2 {x) 


= r(x) 


(ai (x) cn (x) w 3i 2 


/ \ / \ \ k 7 ,i fx) 
- c 2 fx) c 25 fx) u 6>5 ) 

^7,1 


(20) 


k m {x) 


= r(x) 


(61 (x) cn (x) cj 5i6 


/ \ / \ \ ki 7{X) 
C 2 fx) C 52 {X)U 2 ,3) 


(21) 


k m (x) 


= r(x) 


(61 (x) cn (x) cj 6i5 


( \ ( \ -,k 7 \{x) 
C2 \x) c 52 fx)u; 3i2 ) 

^7,1 


(22) 



where we have chosen the entry k\^ 7 {x) as an overall normalization. We also have made use 
of the identities Ci2(x)c2i(x) = Ci 3 (x)c3i(x) = C25(x)cs2(x) in order to define the common 
function, 

= ci (x) cn (x) - c 2 (x) c 4 (x) 

c 2 u (x) b x (x) ai (x) - c\ (x) C12 (x) c 2 i (x) 



At this point we turn our attention to the diagonal B[i,i] blocks. Each of them involves 
two distinct equations and after some cumbersome manipulations we find the relation, 

wp,a k a A x ) = v<x,p kp >a (x), for a ^ (24) 
provided that the following constraints between the coefficients uj a ^ are satisfied, 

W a ,j9 = Wp >a UJ/3>, a >, for^«,«' (25) 

The last off-diagonal elements we need to determine are the entries k^p{x) for (3 > 4 and 
k a ^{x) for a < 4. It turns out that they can be fixed by using the blocks B[l, 17], B[9, 33], 
5(2,20], B[3, 12] and 5(8,35]. The final results are 

fci,4(a0 = — k 2 ,a(x), k 2A (x) = -^k li5 (x), k 3A {x) = — k lfi (x) (26) 

^2,3 ^1,5 w l,6 

k 4 ,5i x ) = — k2j(x), h, 6 {x) = ^-k 3J (x), k AJ {x) = — Afc^a?) (27) 

^2,7 ^3,7 ^5,6 

We now reached a point in which all the above considerations can be collected together. 
Substituting the determined off-diagonal elements back to the reflection equation (2) we find 
out that the parameter u)2j has to satisfy the following polynomial equation, 

^2j(^2,7 - ^1,6<T 5 ) (^2,7 + ^1,6<T 5 ) (W 2 ,7 + ^1,6<T 8 ) = (28) 

The values 002,7 — 0, —oo\fiq~ % lead us to solutions whose non-null elements are only the 
diagonal k a ^ a {x) and the anti-diagonal fc Qja / entries. This branch gives origin to two distinct 
classes of solution depending on whether cj 1i7 = or cj 1)7 7^ 0. By setting ou 1 j = we get the 
two diagonal solutions found previously in ref. [15], which in current notation reads, 

K^\x) = h (29) 

K^\x) = Bmg(l,l,x^^,x^^,x^^,x 2 ,x 2 ) (30) 

q + ex q + ex q + ex 

where e = ±1 is a discrete parameter. 
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On the other hand, for u±j ^ we obtain the following novel non-diagonal solution 



/ 1.(3) 



K^\x) 



kf}{x) 



\ 



k 



(3) 
2,2 



1.(3) 
ft 2,6 



1.(3) 
ft 3,3 



1.(3) 
ft 5,3 



4,4 



1.(3) 
ft 3,5 



.(3) 
"5,5 



1.(3) 
H,2 



U(3) 
v 6,6 



where the diagonal entries are given by, 



K 7,7 J 



1.(3) 



1.(3) 
K44 



,(3) 
v 2,2 



1 g 2 



,(3) 
b 3,3 



2 2 



7 q 2 — I x 2 



1,(3) _ . (3) _ , (3) 
ft 5,5 _ ^6,6 _ ^7,7 

while the off-diagonal elements are, 



f 2x 



—jk[j(x) 



UJlj x z 



k {3) 

^2,6 — 



1.(3) . 
ft 3,5 — 



k 



(3) 
5,3 



^1,7 
1 

^2,6 

^2,6 
_ , ,2 



2g 



2g 



q 2 - 1 



1.(3) 

^6,2 — 



1.(3) 
ft 7,l 



1 / 2g ' 

^1,7^2,6 V? 2 _ 

2 



2g 



(31) 



(32) 
(33) 
(34) 



(35) 
(36) 
(37) 
(38) 
(39) 



^1,7 W ~ 1 , 

The remaining branch u 2 j = ec^i j6 g~ 5 is the most complicated one since it leads us to K- 
matrices with all non-null elements. Denoting the entries of this complete matrix by k^ a ^{x) 
we find, after elaborated algebraic manipulations, that their expressions can be defined as 
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follows. The diagonal entries are, 



q 9 (l + q) (q 8 + ex)(l + q 4 ) ] J4 ) 



(q 3 + 1) (x 2 - 1) (q 5 - ex) 



kf> (x) 



(40) 



= kf\ (x) - J#> (2, 2) q 5 (9 _ +1) M 4 j (*) 



g J — ex 



(41) 



fcgto = fcg(x)+^>(2,2) gl ° s (g6 - 1) JbW 



g 5 — ex 



x 



(42) 



fc W^ _ We) f2 2 x / - * 2 ) _ g 9 (q 9 + l)(eqx + l) \ (4) 

fc 4 , 4 (x)-.A ^^( v(g3 + 1)(;r2 _ 1) (;r 2_ 1)(g 5_ e;r) "jVW 



(43) 



fcg(x)=x a fcS5(x)+^ (2,2) 



,2 7.(4) 



^-1)^(4) 



g° — ex 



ex/c} 7 (x) 



(44) 



x 2 Aff(a:) + .^(2,2) ^ + 



g 5 — ex 



(x) 



(45) 



= x^ (x) 



2/, (4) 



(46) 



The structure of the off-diagonal entries shall be described in terms of the elements of each 
row separately. For the first row we find, 

,9 



l(4),V 

&i 4 ^ (x 

1,(4)/ 
«, 15 ^X 

^16^ 



= -j^(l,2) 



5 



(x) , 



g J — ex 

= (2, 1) (*) , 

q b — ex 



= -^ (e) (0,l) 



5 



g J — ex 

_^)(i >0 )^_ 
g° — ex 



q° — ex 
k$ (x) , 



(47) 
(48) 
(49) 
(50) 
(51) 



The second row is, 



k 2 \ (x) 

The third row is, 

k$(x) 

The fourth row is, 



^(3,2)-^ (52) 
(3, 1) -j^— A$ (x) , (53) 

^10 

g u — ex 



q 2 (1 + ^(e) (2> 1} fc W (x) j (54) 

„10 

-J*> (1, 1) fc$ (x) , (55) 



g 5 -ei 

^5 



-^0(2,0)^ (56) 
-^)(1,0)^— exk^ix). (57) 



-J*>(2,3)^ *$(*), (58) 

q° — ex 

„16 

^)(l,3)-^-^(^), (59) 

-V? 2 (l + ? 2 )^ (e) (1, 2) -r— (*) , (60) 

^)(0,2)-^£$(*), (61) 

-J*> (1, 1) -y^— ex*$ (x) , (62) 
g° — ex 

J*> (0, 1) exkf} (x) . (63) 
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k { £{x) 

^43 ( x ) 

k{$ (x) 
k ( S(x) 

The fifth row is, 

1,(4)/ 

1,(4)/ 
^52 l x 

^54 ( x 

t,(4)/ 
^56 \ x 

1,(4)/ 

The sixth row is, 



-VV 2 (l + <7 2 )^ (£) (3, 3) A$ (*) , 

v q° — ex 

y/q*(l + q*)jM (2, 3) -J^— fc$ (x) , 
v q° — ex 

Vg 2 a + ^ (£) (3,2)-^-A;g(x), 

-y/q2(l + q2)jM (1,2) - T ^— ca *$ (*) , 
v q b — ex 

-VV(l + ? 2 )^ (£) (2, 1) cxM? (*) , 

vV(i + g 2 )^ (£) (i, i) ^X- e*A$ (*) • 

24 

-•^ (4,3)^^*5S(x), 
q b — ex 

23 

(3, 3) fc$ (*) , 

q b — ex 

„18 

V(l + g 2 )g 2 ^ {e) (3, 2) - T ^!_ ca *$ (x) , 

( 3 ,l)^^ e ^(x), 
q b — ex 

J** (2,1)^— e xfcJ§(x). 
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h$(x) = .^(3,4)^— M?t>), (76) 
k$(x) = -&<)(2,4)-^k§(x), (77) 

„18 

*$?(*) = ^ (e) (3,3)^ exk§(x), (78) 

q° — ex 

*£>(*) = + g 2 ) tfj™ (2, 3) exk$ (x) , (79) 

= ^ (e) (l,3)^^6a;A;S 4 )(x), (80) 
= -P^(l,2)^—exk{ 4 }(x), (81) 

The seventh row is, 

„26 

fcff(x) = J*>(4,4)-^ *$(*), (82) 

„22 

fcg)(x) = ^)(3,4)-^-— (83) 

22 

fcff(x) = J*)&Z)^—e X k§{x), (84) 
h$(x) = -VPa + g 2 )^ (e) (3, 3) - T ^!_ esfc$ (*) , (85) 

„19 

= (2, 3) exfc^(x), (86) 

„19 

= J*> ( 3 , 2) ^r^— (87) 
where the auxiliary function jF( e ) (n, to) is defined by, 

on+m / „3 I i \"+m+l 

77(6) / N = f 19 ++J /ggN 

(e + g) (e + g 4 ) (e + g 8 ) ^ 2 uj\% 

Note that this latter solution contains an additional discrete parameter besides the two 
continuous uj\^ and o; li3 variables. We also observe that a striking feature of our results is that 
the non-diagonal solutions do not reduce to any of the admissible diagonal matrices K^>{x) 

11 



and K^(x). This means that the four possibilities discussed above are indeed the distinct 
types of if-matrices of the U q [G 2 } model. The latter feature should be contrasted with the 
findings for the non-exceptional vertex models [11] in which reductions of general non-diagonal 
matrices guide us always to some particular diagonal solution. It remains to be seen whether 
such peculiarity is special to the U q [G<2\ symmetry or is valid for all exceptional vertex model 
as well. 

Finally, equipped with the four reflection matrices K^(x), an integrable model with open 
boundary condition can be obtained through the double-row transfer matrix formulated by 
Sklyanin [3], namely 



where T(x) = R a L(x) ■ ■ ■ R a i(x) is the standard monodromy matrix of the corresponding closed 
chain with L sites. The matrix (x) is automatically determined [5] from K ( m ) (x) with the 
help of the crossing property (7), 



In this paper we have been able to classify the possible families of reflexion matrices associated 
to the fundamental U q [G 2 ] vertex model. We find that there exists four different classes of K- 
matrices. Two of them are the diagonal solutions studied before [15] and the remaining ones are 
new non-diagonal if-matrices with continuous free parameters. The first type of non-diagonal 
solution is made by combining only the diagonal and the anti-diagonal elements while the 
second one consists of matrix whose all entries are non-null and has an extra discrete parameter. 
Interesting enough, these non-diagonal solutions do not possess reduction to particular diagonal 
matrices for any values of the free parameters at disposal. We stress that this peculiarity does 
not occur in the ^-matrices of the vertex models based on the non- exceptional Lie algebras 



t^ m \x) = Ti a K + (x)T(x) K (x)T-\l/x) for I, m = 1, • • • , 4 



(89) 



K^\x) = \K^ m \q 12 /x)] V l V 



(90) 



4 Conclusions 



[ii]. 
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We hope that our results will prompt further lines of investigations. One possibility is 
to explore the Weyl representation of the U q [G 2 ] -R-matrix to establish, at least for diagonal 
boundaries, the diagonalization of the transfer-matrix (89) from a first principle framework 
such as the quantum inverse scattering method [16]. This will give us an opportunity to check 
certain assumptions, usually denominated "doubling hypotheses" [17], used to determine the 
corresponding Bethe ansatz solution [15]. . 
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Appendix A: The Boltzmann weights 



In this Appendix we are going to describe all the Boltzmann weights appearing in Eq.(5). In 
order to do that we find convenient to introduce an auxiliary function ip(a,(3) defined by, 



<p(a,P) 



[f ] - a + 4<5 4i/3 



for a < (3, a < (3' 
for a < [3, a > f3' 
for a > j3, a' > f3 
for a > (3, a' < (3 



(A.I) 



[f]-[3' + 45 a , A 
/3-[f]-4<f a , 4 
( a' - [f ] - 4 V 

where [|] denotes the integer part of |. 

We begin by listing the weights Ai(a,j3), A^{a) and A[k] as well as Bi(a,{3), 82(a) and 
B[k]. They are given by, 



Ai{a,(3) 



a 7 (x) for (p(a, (3) = —1 

for <p(a, (3) = ,A 2 {a) = 
a 6 (x) for (p(a, (3) — 1, 2 
a 2 (x), A[ 2 ] = a 3 (x), A[ 3 ] = a 5 (x) 



04(2;) for (p(a, 4) = 5 
a 8 (x) for (p(a, 4) 7^ 5 



, (A.2) 



(A.3) 
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B 



[i] 



&7(x) for (p(a,(3) = — 1 

61 (x) for <p(a,P) = ,i3 2 (a) = 

6 6 (x) for (p(a,P) = 1,2 

6 2 (x), #[2] = 63 (») , S[3] = 



64(0;) for (/?(«, 4) = 5 
6s 0*0 for </?(«, 4) 7^ 5 



while the weight C(a, (5) is, 



c 2 (a;) for (p(a, (3) — 
04(0;) for </?(ai, /3) 7^ 



C(a,/3) = 

The functions a a (x) and 6 a (^) entering in the above expressions are, 



ai(u 
a 2 (u 

a 3 (u 
a 5 (u 

a^x 

a & (x 

ar(x 

a$(x 



where 



(g 2 - l) (g 8 - x) (g 12 -x), b± (x) = xcn (x) , 
-q 7 y/l + <T 2 {q 2 ~ l) (q 12 -x)(x-1), b 2 (x) = fax), 
q 4 \Jl + q~ 2 (q 2 - l) (g 12 - xj (x - 1) , b 3 (x) = xqa 3 (x), 
q^l + q- 2 (q 2 - l) (<? 12 - x) (x - 1) , 6 5 (x) = fax), 

(g 4 - l) (g 12 - x) fi(x), b 4 (x) = ^r^- a 4 (x), 

Ji\ x ) 

(q 2 - l) (g 12 - x) g^x), b 6 (x) = ^^a 6 (x), 
(q 2 - l) (Y 2 - x) g 2 (x), b 7 (x) = 
(<? 4 - l) {q 12 ~ x) f 2 (x), b 8 (x) = J^*s(xl 



h(x) = g 2 (l-g 2 + g 4 )-x, 

f 2 (x) = g 6 - (l - g 2 + g 4 )x, 

^(x) = g 2 (l + g 2 + g 6 )-(l + g 2 + g 4 ):r, 

g 2 (x) = q\l + qi + q*)-(l + q 4 + q«) x . 
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Next, the weights c a (x) are, 

ci (x) = (g 2 - x) (q 8 - x) (q 12 - x) , 

c 2 (x) = - q ( q 8 - x )(q 12 - x )( x -l), 

c 3 (x) = -q 2 {q 6 - x)(q 12 - x)(x - 1), 

c 4 (x) = - X )( q » - x )(x - 1) 

while the c Q/ 3(x) weights are given by 

-q 4 (q 4 - x) (q 10 - x) (x - 1) , (a = (3, a ^ 4) 
( ? 6 -x)[x 2 + (l + g 4 )[( ? 4 -l) 2 (? 2 -l) 2 - g 6 ]x + g 14 ], (a = /3 = 
CQ/3 X ~ ' (-l)^(g 2 - l)g 16 +^(:r - l)^ K/3) , (a </?, a') 

where the auxiliary index a is defined as, 

1 = 0, 2 = 3, 3 = 6, 4 = 8, 5 = 10, 6 = 13, 7 = 16 
and the functions ^>( Qj/ 3) are 

^"o = g 4 - a; 

^ 5 = ^ 3 = ^-4 = ^(l + g 2 )[g 6 (l-g 2 + g 4 )-:r] 
= ^ 3 = = 1(1 + g 2 ) [g 10 - (1 - g 2 + g 4 )x] 

Tx = ^ 2 = l[g 8 (l + g 4 + g 6 )-(l + g 2 + g 4 )x] 
= ^ 2 = l[g 10 (l + g 2 + g 4 )-(l + g 2 + g 6 )a;] 
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Finally, the remaining weights c Qja /(x) are given by 



c 17 (x) 


= (q 2 - 


l)(g 8 -l)[g 8 (l + g 4 )- 


(q 8 + l)x] 




(A.30) 




= (q 2 - 


1)[? M (? 10 -1)-?V 


- l)(q 4 + l)x - (q 6 - 


l)x 2 ] 


(A.31) 




= (a 2 - 




- l)x - (q 12 - l)x 2 } 




(A.32) 




= (s 2 - 


i)[?V 2 -i)-?V- 


- l)x - (q 4 - l)x 2 ]x 




(A.33) 




= (V 2 - 


l)[q U (q 6 -l)+q 8 (q 2 - 


-l)(q 4 + l)x-(q 10 - 


l)x 2 ]x 


(A.34) 


c 71 (x) 


= (? 2 - 


i)(? 8 -i)[?V + i)- 


(q 4 + l)x]x 2 




(A.35) 
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